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Abstract 

We examine the potential of the COMPASS experiment at CERN to study color 
transparency via exclusive coherent vector meson production in hard muon-nucleus 
scattering. It is demonstrated that COMPASS has high sensitivity to test this important 
prediction of perturbative QCD. 
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1 Introduction 



One of predictions of perturbative QCD (pQCD) is that at high Q 2 the longitudinally 
polarized virtual photons 7^ fluctuate into hadronic components, e.g. qq pairs, whose 
transverse size b = \ f± q — f±q | decreases with Q 2 , b 2 oc (1/Q 2 ). At large Q 2 the values 
of b are significantly smaller than the size of the nucleon. Color transparency (CT) is 
a phenomenon of pQCD, whose characteristic feature is that small color-singlet objects 
interact with hadrons with small cross sections [1, 2, 3]. Cross section for the interaction 
of such small object, or small size configuration (SSC), with a hadron target has been 
calculated in pQCD using the factorization theorem [4, 3, 5], as: 

Vqq, N = y b 2 a s {Q 2 ) x g(x, Q 2 ) , (1) 

where g(x, Q 2 ) is the gluon distribution function in the nucleon. 

Observing CT in particular kinematics would prove experimentally the applicability 
of the pQCD factorization theorem for those kinematics. Such a proof provides impor- 
tant complementary support to a class of spin physics experiments, for example related to 
measurements of generalized parton distributions at COMPASS [6]. 

The prerequisite for observing CT is to select a sample containing SSC's via a hard 
process (i.e. with large Q 2 , high p t , or large produced mass). For hard exclusive p° lep- 
toproduction, in addition to large Q 2 , selection of the longitudinally polarized mesons is 
required. To probe the SSC-nucleon cross sections one should study absorption of the SSC 
propagating through nuclear matter. In order to clearly observe CT it is necessary that 
the SSC lives long enough to traverse distances larger than the size of the target nucleus, 
and that the SSC stays small while propagating through the nucleus. These requirements 
are characterized in terms of the coherence length l c (the distance traversed by the qq fluc- 
tuation of the virtual photon in the target nucleon/nucleus system) and by the formation 
length (the distance in the target system after a scattering needed for the qq fluctuation to 
develop into a hadron). 

We consider here only exclusive coherent vector meson V production at small \t\. We 
analysed the coherent and incoherent channels in Refs [7, 8]. For near-forward (t « 0) 
coherent production and complete CT, one expects that the cross section for the nuclear 
target is related to that for the nucleon target, provided that |i m i n <R\> /3| <C 1, by: 

^ (7 M -> VA) U = A 2 ^*N - VN) U . (2) 

Usually in experiments the t-integrated coherent cross section is measured, for which 
CT predicts an approximate A 4 / 3 dependence. Without CT, if a hadron propagates in 
nuclear matter, the expected ^-dependence is weaker; A 2 / 3 for the cross sections comparable 
to the pion-nucleon cross section. 

In searches for CT a commonly used quantity, measured in experiments, is the nuclear 
transparency 
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which is the ratio of the cross section per nucleon for a selected process on a nucleus A to 
the corresponding cross section on a free nucleon. 

We concentrate on exclusive coherent vector meson production (EVMP), which could 
be studied at the COMPASS experiment [9, 10]. Kinematic variables used in this paper are 
listed in Table 1. 
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Table 1: Kinematic variables used in the text. 



k, k 




four-niomcntum of the incident and scattered muon, 


p, V 




four-momentum of the target nucleon and vector meson V, 


i = 


k-k' 


four-momentum of the virtual photon, 


Q 2 


= -1 2 


virtuality of the photon 7*, 


v = 


(p ■ q)m, 


energy of the virtual photon in the laboratory system, 
Mp is the proton mass, 


x — 


Q 2 /(2M pV ) 
{q-vf 


Bjorken scaling variable, 


t = 


squared four-momentum transfer to the target, 


P 2 




transverse-momentum squared of the vector meson, 


my 




invariant mass of the vector meson V, 


Ml 


= {p + q-v) 2 
(M| - Mp 2 )/W 2 


missing-mass squared of the undetected recoiled system, 


I = 


inelasticity 



In order to study CT for the exclusive production of mesons composed of light quarks 
(like p) one should select large Q 2 , moderately small x and longitudinally polarized mesons. 

The t-dependence of the cross section for exclusive meson production on the nucleon 
is approximately exponential, e B ^~ tmlTi \ where i m i n « —M 2 x 2 (l + rrty/Q 2 ) 2 and my is 
the mass of the produced meson. For hard EVMP, if SSC's are important, the slope B 
significantly decreases with increasing Q 2 , approaching a universal value (4.5 — 5 GeV~ 2 ) 
related to the nucleon radius. This is so because the transverse interquark distance in the 
SSC decreases with increasing Q 2 . 

For EVMP on a nucleus A the t-dependence of the cross section is approximately 
reproduced by a sum of two exponential functions. The peak at the lowest \t\ values, 
with the slope proportional to the nucleus radius squared, < R\ >, is mostly due to the 
coherent production, whereas at somewhat larger \t\ the incoherent production on quasi- 
free nucleons dominates and the slope is almost equal to that for the production on free 
nucleons. Applying cuts on t allows to select successfully samples of events, which are 
strongly dominated by either coherent or incoherent production, as discussed in Refs. [11, 
12, 13]. 

Strong recent evidence for CT comes from Fermilab E791 experiment on the ^-de- 
pendence of coherent diffractive dissociation of pions into two high-p t jets [13]. Also the 
E691 experiment results on ^-dependence of coherent J/tp photoproduction [12] and the 
NMC measurements of ^.-dependence of coherent and quasielastic J/0 muoproduction arc 
consistent with CT, but don't demonstrate it unambiguously. 

2 Experimental method 

We proposed [7] to study CT via exclusive coherent p meson production pA — > p pA on 
C and Pb nuclei, as well as the exclusive production of other vector mesons. The preferable 
decay modes are those into the charged particles such as p° — > 7r + 7r~. 

We propose to complement the initial setup of the COMPASS [9, 10] for the muon- 
beam run with the polarized target by adding the nuclear targets downstream of the polar- 
ized target. 

The high-intensity high-energy incident muon-beam will impinge on the polarized 
target and a downstream thin nuclear target. The momenta of the scattered muon and of 
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the produced charged particles will be reconstructed in the two magnetic spectrometers, 
using the magnets SMI and SM2, instrumented with various tracking detectors. In order 
to select exclusive events, we rely only on the kinematics of the scattered muon and the 
produced meson. 

The trigger will use the information from hodoscopes registering the scattered muons, 
and from calorimeters registering deposits of energy of the particles in the final state. 

The discrimination of non-exclusive events will be done by applying cuts on the in- 
elasticity I (for the definition see Table 1). In Fig. 1 the inelasticity distribution is shown 
for the SMC p° sample [14] for the events with the invariant mass in the central part of 
the p° peak. For the inelasticity distribution the peak at I = is the signal of exclusive p° 
production. Non-exclusive events, where in addition to detected fast hadrons, slow unde- 
tected hadrons were produced, appear at / > 0. With larger statistic in COMPASS it will 
be possible to apply more tight inelasticity cuts in order to further reduce the background. 

The selections of coherent or incoherent production will be done on a statistical basis, 
using the t-distribution; at the lowest \t\ values coherent events predominate. 

Separation of the p° samples with the enhanced content of longitudinally or trans- 
versely polarized mesons will be done by applying cuts on the measured angular distributions 
of pions from the decays of the parent p°. 

The minimal covered Q 2 -range is expected to be 0.05 < Q 2 < 10 GeV 2 . For the med- 
ium and large Q 2 values (Q 2 > 2 GeV 2 ) the range 0.006 < x < 0.1 will be covered with 
good acceptance. 

The basic observable for each process studied will be the ratio of the nuclear trans- 
parencies for lead and carbon, R T = T^/Tq = (<Tpb/^4pb)/(fcMc)- 

3 Simulation of exclusive p° events 

We present details of the simulation of exclusive coherent p° production (pA — ► p p°A) 
in the COMPASS experiment with the carbon and lead targets. The simulations were done 
with a dedicated fast Monte Carlo program which generates deep inelastic exclusive p° 
events with subsequent decay p° — > ir + ir~. At this stage there was no attempt to include 
any background in the event generators. Details of the used parameterization of the cross 
section for the production of p° on the free nucleon, pN — ► p p°N, with the subsequent 
decay p° — > 7r + 7r~ were presented in Ref. [7]. We relate the differential cross sections for 
the free nucleon to these for coherent production on the nucleus A via 

( da A h ) _ A 2 r <Rl>t/3 ( d<r N ) (4) 

Here < R\ > is the mean squared radius of the nucleus, A e ^ co j 1 takes account of nuclear 
screening for the coherent process, and i = L or T designates the polarization of p°. We 
used the approximation t — t m i n ~ — p 2 . 

We generated the cross sections for two models. For the complete color transparency 
model (CT model) we used A eS coh = A. In another model we assumed a substantial nuclear 
absorption (NA model) and used A eS coh = A ' 75 . 

Details of the simulation of the apparatus effects were given in [7]. Here we discuss 
only the trigger simulation. To simulate the trigger acceptance a trajectory of the scattered 
muon behind the second magnet was calculated, and the hits in the muon hodoscopes 
were checked. We consider two different trigger acceptances, the Medium Q 2 range Trigger 
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(MT) and the Full Q 2 range Trigger (FT). The Q 2 -dependence of the trigger efficiency is 
denoted e tr . For the MT trigger at 190 GeV beam, e tr decreases quickly with Q 2 from 0.7 
at Q 2 = 2 GeV 2 to about 0.1 at Q 2 = 10 GeV 2 For the FT trigger the Q 2 -dependence is 
weaker and the trigger acceptance is higher; it is always bigger than 0.5 for Q 2 < 70 GeV 2 
at 190 GeV beam energy. 



4 Results on exclusive p° production 

We considered 190 GeV muon beam. The simulations were done independently for 
the carbon (A = 12) and lead targets (^4 = 207), and for two triggers (MT and FT). For 
each target and each trigger we assumed two different models describing the nuclear effects 
for exclusive p° production: CT model and NA model (cf. Section 3). 

The kinematic range considered was the following: 

2 < Q 2 < 80 GeV 2 , (5) 

35 < v < 170 GeV . (6) 

The upper cut on v was chosen to eliminate the kinematic region where the amount of 
radiative events is large, whereas the lower one to eliminate the region where the acceptance 
for pions from p° decay is low. 

The total efficiency e tot to observe exclusive p° events results from: the acceptance of 
the trigger (e tr ), acceptance to detect the pions (e had ), efficiency for tracks reconstruction 
(e rec ), cut on the invariant mass of two pions, 0.62 < < 0.92 GeV, used for the selection 
of the samples (e mass ), and efficiency for an event to survive the secondary interactions (e sec ). 
The contributions of all these effects to e tot are similar to those presented in [15] for the 
polarized target, except for the effects of the secondary interactions. The approximate value 
of e sec is equal to 0.87 for the carbon target and 0.94 for the lead target. The total efficiency 
e tot is about 0.48 for the FT trigger and about 0.30 for the MT trigger. 

The total expected cross section for exclusive p° production on the nucleon is 



r^max f'Q max _ 

°Tn^ P °n = J J %n^ p °n(Q > ") dC ? dv > ( 7 ) 

min ^min 

where the kinematic range was defined before. The value of ^J^v pojv is 283 pb. For 
nuclear targets the corresponding values depend on A and on the assumed model for nuclear 
absorption. 

The expected muon beam intensity will be about 10 8 /s during spills of length of about 
2 s, which will repeat every 14.4 s. With the proposed thin nuclear targets, each of 17.6 
g/cm 2 , the luminosity will be C = 12.6 pb -1 • day -1 . The estimates of the numbers of 
accepted events were done for a period of data taking of 150 days (1 year), divided equally 
between two targets. An overall SPS and COMPASS efficiency of 25% was assumed. The 
numbers of accepted events for the carbon and lead targets, assuming the two models for 
the nuclear absorption mentioned earlier, are given in Table 2. 

The distributions of accepted exclusive events as a function of x, Q 2 and v are similar 
to those presented in [15]. 

In Fig. 2 we present the p 2 distributions for both targets. We observe clear coherent 
peaks at small p 2 (< 0.05 GeV 2 ) and less steep distributions for the incoherent events at 
larger p\. The arrows at the top histograms indicate the cut p\ > 0.1 GeV 2 , used to select 
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Table 2: Numbers of accepted events for two considered models of the nuclear absorption. 



model 


N c 




CT 


70 000 


200 000 


NA 


28 000 


20 000 



the incoherent samples. For the middle and bottom histograms the arrows indicate the cut 
pi < 0.02 GeV 2 , used to select the samples which are dominated by coherent events — the 
so called coherent samples. For the samples defined by the latter cut the fraction of the 
incoherent events is at the level of up to 10%, depending on the nucleus and on the model 
for nuclear absorption. 

The effect of the kinematical smearing on p 2 may be seen by comparing the distri- 
butions for the generated events (middle row) to the ones for measured events (bottom 
row) where the acceptance and smearing were included. The smearing of pi increases with 
increasing p 2 ; it is about 0.006 GeV 2 for the coherent samples {pi < 0.02 GeV 2 ) and about 
0.03 GeV 2 for the incoherent samples (pi > 0.1 GeV 2 ). 

The analysis of the p° decay distributions allows us to study spin-dependent proper- 
ties of the production process [16], in particular the polarization of p°. The p decay angular 
distribution W(cos 9, 4>) is studied in the s-channel helicity frame. The p° direction in the 
virtual photon-nucleon centre-of-mass system is taken as the quantization axis. The angle 
9 is the polar angle and </> the azimuthal angle of the 7r + in the p° rest frame. The cos 9 
distributions for pions from p° decays are shown in Fig. 3 for the lead target. The distribu- 
tions for longitudinally (dashed lines) and transversely (dotted lines) polarized parent p°'s 
are markedly different. Their sum is also indicated. 

As Q 2 increases, the approach to the CT limit is expected to be different for EVMP 
by longitudinally polarized virtual photons from that by transversely polarized photons. To 
increase the sensitivity of the search for CT, we study ^-dependence of the cross sections for 
samples with different p° polarizations, which will be selected by cuts on cos 9. Such cuts 
will allow us to select the samples with enhanced contributions of the events initiated by the 
virtual photons of a desired polarization. Studies of the samples with different polarizations 
of the virtual photons are necessary for the unambiguous demonstration of CT. 

Another aspect which is important for CT studies, is the covered range of the coher- 
ence length l c (cf. Section 1). We showed in Ref. [7] the correlation of l c vs. Q 2 for a sample 
of accepted events. At small l c values the effects of initial and final state interactions in the 
nuclei mimic CT [17]. It is important to disentangle effects due to CT from those caused 
by the modified absorption at small l c values. 

For a limited statistics, we may use the combined data at l c values exceeding the 
sizes of the target nuclei. The radius of the carbon nucleus is < r 2 ^ >V 2 = 2.5 fm and 
that of the lead nucleus is < rp b > 1 / 2 = 5.5 fm. Therefore, one may use the selection 
l c > ~ 2- < rp b > l / 2 = 11 fm. About a half of events survive such a cut on l c . These 
events cover the range of Q 2 < 6 GeV 2 , which is expected to be sufficient to observe CT. 

The estimated values and statistical precision of i? T , the ratio of the nuclear trans- 
parencies for lead and carbon, are presented for different Q 2 bins in Fig. 4 for pi < 0.02GeV 2 . 
The Q 2 bins are specified in Table 3. Figure 4 comprises predictions for two models, CT 
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Table 3: Q 2 bins used for the determination of R T . 



Kill VIllTYlhDT 1 

J_>111 IlUIilUt:! 


D 2 Kin 
\cg Dill 


<r n 2 -> 






[GeV 2 ] 


[GeV 2 ] 




1 


2-3 


2.4 


0.016 


2 


3-4 


3.4 


0.022 


3 


4-6 


4.8 


0.031 


4 


6-9 


7.2 


0.048 


5 


9-12 


10.2 


0.072 


6 


12-20 


14.8 


0.11 



and NA, and for 6 different samples of accepted events for each model. For each sample 
a set of "measurements" in different Q 2 bins is shown. Sets A and B were obtained us- 
ing the standard selections for the MT and FT triggers, respectively. Four remaining sets 
were obtained for the FT trigger with additional selections: C with | cos# |< 0.4, D with 
| cos6> |> 0.7, E with | cos 6* |< 0.4 and l c > 11 fm, and F with | cos# |> 0.7 and l c > 11 fm. 
Note that for sets E and F only three lower Q 2 bins appear. One expects large differences 
in R T for the two considered models. For coherent samples R T f» 5 for CT model and « 1 
for NA model. At Q 2 ~ 5 GeV 2 the precision of the measurement of i? T for coherent events 
will be better than 17%, even for the restricted samples E and F, thus allowing excellent 
discrimination between the two models for the nuclear absorption. 

5 Comparison with previous experiments 

We compare COMPASS capabilities to demonstrate CT to those of previous experi- 
ments in which exclusive coherent p° leptoproduction on nuclear targets was studied. For 
the comparison we must also select experiments which covered Q 2 range extending to large 
values, bigger than 2 GeV 2 . The only experiment which published such high Q 2 exclusive 
coherent data is the NMC experiment [18]. 

NMC has published data on coherent exclusive p° production on 2 H, C and Ca targets. 
The muon beam energy was 200 GeV and the data cover the ranges 2 < Q 2 < 25 GeV 2 
and 1 < l c < 30 fm. Due to the moderate statistics of the data the Q 2 -dependence of 
nuclear absorption was not obtained at sufficiently large or at fixed l c , which is necessary 
to unambigously demonstrate CT. 

The COMPASS data at medium and large Q 2 (FT trigger) will cover the kinematic 
range similar to that of the NMC data. The expected statistics for the carbon target will 
be about 2 orders of magnitude higher than that of the NMC data for the same target. 
This increase is in particular due to higher beam intensity and larger acceptance in the 
COMPASS experiment. In addition COMPASS will use also lead target and will extend 
measurements to small Q 2 . Due to large statistics splitting of COMPASS data in several Q 2 
and l c bins as well as the selection of events with longitudinal or transverese p° polarization 
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will be possible. 



6 Conclusions 

We have simulated and analysed exclusive p° muoproduction at the COMPASS exper- 
iment using thin nuclear targets of carbon and lead. For muon beam energy of 190 GeV 
and a trigger for medium and large Q 2 , the covered kinematic range is 2 < Q 2 < 20 GeV 2 
and 35 < v < 170 GeV. 

Good resolutions in Q 2 , l c , t (p 2 ) and cos 9 are feasible. An efficient selection of 
coherent events is possible by applying cuts on p 2 . In order to obtain the samples of events 
initiated with a probability of about 80% by either or 7^, the cuts on the p° decay 
angular distribution of cos 9 will be used. The search for CT could be facilitated by using 
the events with l c values exceeding the sizes of the target nuclei. The fraction of such events 
is substantial and the covered Q 2 range seems sufficient to observe CT. 

We showed high sensitivity of the measured ratio R T of nuclear transparencies for 
lead and carbon for different models of nuclear absorption. Good statistical accuracy of the 
measured i? T may be achieved already for one year of data taking. These measurements, 
taken at different Q 2 intervals, may allow to discriminate between different mechanisms of 
the interaction of the hadronic components of the virtual photon with the nucleus. 

In conclusion, the planned comprehensive studies of exclusive vector meson produc- 
tion on different nuclear targets at the COMPASS experiment are able to unambiguously 
demonstrate CT. 
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Figure 1: Inelasticity distribution. The SMC preliminary results [14] for fiN — > fip°N. 
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Figure 2: Distributions of p\ for the lead (left) and carbon (right) targets. The distributions 
a) and b) at the top are for the generated events and the wide range of pi , whereas those 
in the middle, c) and d), correspond to the range of low p\. The bottom distributions, e) 
and f), are for the accepted events, with the kinematical smearing taken into account. The 
arrows show the cuts described in the text. 
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Figure 3: Distributions of cos 6 for the pions (tt + ) from the decays of p° mesons produced 
on lead, for generated (left) and accepted (right) events. For the latter ones the effects of 
the kinematical smearing were taken into account. The dashed- and dotted-line histograms 
are for the longitudinally and transversely polarized p°'s, respectively, and the solid-line 
histograms are for their sum. 
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Figure 4: Expected statistical precision of the measurement of the ratio R T of GeV 2 nuclear 
transparencies for lead and carbon for pi < 0.02 GeV 2 and for different Q 2 bins, for CT 
(upper band) and for NA (lower band) models. The Q 2 bins are defined in Table 3. Sets 
A, B, C, D, E and F correspond to the conditions specified in the table below the plot. 
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